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ABSTRACT 

The origin of the extragalactic 7 -ray background (EGB) has been debated for 
some time. The EGB comprises the 7 -ray emission from resolved and unresolved 
extragalactic sources, such as blazars, star-forming galaxies and radio galaxies, 
as well as radiation from truly diffuse processes. This letter focuses on the blazar 
source class, the most numerous detected population, and presents an updated 
luminosity function and spectral energy distribution model consistent with the 
blazar observations performed by the Fermi Large Area Telescope (LAT). We 
show that blazars account for % of the EGB photons (>0.1GeV), and 

that Fermi -LAT has already resolved ~70 % of this contribution. Blazars, and in 
particular low-luminosity hard-spectrum nearby sources like BL Lacs, are respon¬ 
sible for most of the EGB emission above 100 GeV. We find that the extragalactic 
background light, which attenuates blazars 7 high-energy emission, is responsible 
for the high-energy cut-off observed in the EGB spectrum. Finally, we show that 
blazars, star-forming galaxies and radio galaxies can naturally account for the 
amplitude and spectral shape of the background in the 0.1-820 GeV range, leav¬ 
ing only modest room for other contributions. This allows us to set competitive 
constraints on the dark-matter annihilation cross section. 

Subject headings: cosmology: observations - dark matter - gamma rays: diffuse 
background - galaxies: jets galaxies: active -surveys 


Introduction 


The Large Area Telescope (LAT, Atwood et ah 200 9) on Fermi has recently allowed a 
broadband, accurate, measurement of the extragalactic 7 -ray background (EGB), the inte¬ 
grated emission of all resolved and unresolved extragalactic GeV sources, ch aracter izing its 


intensity over almost 4 decades in energy between 0.1 GeV and 820 GeV(Ackermann et al. 


2014bl . hereafter AC14). At these energies, the EGB spectrum is found compatible with a 


power law with a photon index of 2.32(±0.02) that is exponentially cut off at 279(±52) GeV 
(AC14). Such cut-off, observed for the first time, may be caused by the extragalactic 
background light (EBL; iGould fe Schrederl Il966t IStecker et al.l Il992l) . Yet, this observa¬ 
tion alone is not sufficient to identify which process or source population is responsible 
for the EGB. Specifically, the EGB may encompass the si gnatu res of proces ses g enerating 


a tru ly diffuse background, like intergalactic shocks (e.g. Loeb &; Waxma/n 2000 : Miniatil 


e.g. 


20021 ), -v-ray emission induced by ultra-high energy cosmic rays in inter galactic space (se e 


Bhattachariee &; Siglll2000h . and dark matter (DM) annihilation (e.g. lUllio et al.l 120021) . 
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Besides truly diffuse processes (Alilers & Salvadol 20 


1 ), unresolved point-like sources 


might be responsible for a substantial part of the EGB (jDermer 2007). At high Galactic 
latitudes, Fermz- LAT has detected blazars, star-formin g galaxies, radio galaxies and mil¬ 


lisecond pulsars (lAbdo et al.l l2010bl: iNolan et al.l 120121) . Extensive analyses were recently 


performed to assess the contribution of all these source classes to the EGB. Blazars, which 
constitute the largest po pulation detected by the LAT, were found to contribute ^20 30 % 


of the unresolved EG B (lAbdo et al.l l2010ct iSingal et al.l 1201 2; Harding & Abazai ianl 12012 


Di Mauro et al. 

2014 

), and a larger fraction in some models ( 

Stecker & Venters 

to 

o 

1 — 1 


Star-forming galaxies produce y-rays in cosmi c-ray interactions, with the acceleration o f 
cosmic rays ultimately powered by star formation ([Thompson et al.ll2007l: lLacki et al.ll2014l) . 
Being dimmer but more numerous than blazars, st ar-forming galaxie s might be responsi¬ 


ble for 10-30 % of the 0.1-100 GeV EGB photons ( Fields et al.ll2010t iMakiva et al.l 12010 


Ackermann et al .1120121 ). A similar argument holds for mi saligned AGN (e.g. radio galaxies) , 
which were recently found to produce ~20 % of the EGB (llnouel 201 ll; Di Maur o et al. 2013). 


However, large uncertainties remain for the contributions of the above source classes. 
In this paper we present (in § [2]) improved modeling of the evolution and of the spectral 
energy distributions (SEDs) of blazars, which allows us to quantify their integrated emission. 
We show (in § [3]) that the integrated emission of blazars, star-forming and radio galaxies 
naturally accounts for the amplitude and spectrum of the new EGB measurement over the 
entire 0.1-820 GeV energy range. We then use this information to place constraints on the 
DM annihilation cross section. Throughout this paper, we adopt H 0 =67km s -1 Mpc -1 , 
Da4' = 1-Da = 0.30. 


2. Analysis and Modeling 


The contribution of blazars to the EGB was already estimated in Abdo et al. (2010c) by 
extrapolating the log A-log S below t he Ferrm-LAT detection threshold. T aking advantage 
of recent follow up observations (e.g., IShaw et al.l 120131: lAiello et al.l 120141) . we derive new 


models for the luminosity and redshift evolution of the whole blazar class and of its SED. 
In this sectio n, these mod e ls are c onstrained using blazar data (flu xes, redshifts and p hoton 
indices) from lAbdo et al.l (|2010dl) . the >10 GeV log A-log S' from I Ackerman net al.l (1201311 
and information on the spectral curvature of blazars (lAckerma/nn et ah. |2014d ) to robustly 
estimate the integrated emission of blazars. 
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m 


We rely on the_sample of 403 blazars detected with test statistic^ >50 at |6| >15° 


Abdo et al. (I201 0dh for which a determination of the LAT detection efficiency exists 


f lAbdo et al.l I2010d ). This sample includes 211 BL Lacs, 186 flat-spectrum radio quasars 
(FSRQs) and 6 blazars of uncertain type. All but 109 BL Lacs have a spectroscopic red- 
shift and I Aiello et al.l 1201411 provide redshift constraints for 104 (out of the 109) BL Lacs. 
The fraction of sources with redshift information is ~99% while the incompleteness of the 
sample, due to unassociated sources that might be blazars, is ~10 %. We do not separate 
the two blazar sub-classes (FSRQs, and BL Lacs) because the larger sample allows a better 
determination of the integrated emission from the whole population, and averages over the 
luminosity functions of the two populations in the regime of overlapping luminosities. In 
order to derive the luminosity function (LF), we use the bootstrap Monte Carlo approach 
developed in Aiello et al. (2014) that allows us to treat properly those sources with an im¬ 
precise redshift, providing a robust error estimate. 


We test models of primarily luminosity evolution (PLE), primarily density evolution 
(PDE) and luminosity-dependent density evolution (LDDE). In all these models the sources 
experience both luminosity and density evolution. The LF at redshift z = 0, for sources 
of 0.1-100 GeV rest-frame luminosity L 1 (in erg s^ 1 ), is modeled as a double power law 
multiplied by the photon index distribution: 


$(l 7 , z = o, r) 


dN 

dL^dVdT 


A 

ln(10)L 7 



72 . e -0.5[r-/x(L 7 


l 2 A 2 


( 1 ) 


where /r and a are, respectively, the Gaussian mean and dispersion of the photon index (r) 
distribution. We allow /j, to change with luminosity as: 


/i(L 7 ) = n* + /3 x [log(L 7 ) - 46]. 

(2) 

The evolution is parametrized by the ‘evolutionary factor’ e(z,L 7 ), which 
PDE and LDDE as: 

is applied to the 

<L(L 7 , z, T) = 4>(L 7 , z = 0,r) x e(z, L 7 ), 

(3) 

and to the PLE as: 


<F(L 7 , z,T) — <&(L 7 /e(z, L 7 ), z = 0, T). 

(4) 

For both the PLE and PDE the evolutionary factor is: 


e(z, L 7 ) = (1 + z) kd e z/i , 

(5) 


1 See Abdo et al. ( 2010bh for a definition. 
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with 

k d = k* + t x [log(L 7 ) — 46]. 

For the LDDE0 it is: 


e(z,L 7 ) 


f 1 + z \- pi(L ^ / 1 + z 

\ 1 + z^L^) J \ 1 + z c (L^) 


-P2(i 7 )l 1 




with 


z c {L. y ) = z* • (L 7 /10 48 )“, 
Pi(Lj) = pl + t x (log(L 7 ) - 46), 
p 2 (L 1 ) = p* + 5x (log(L 7 ) - 46). 


( 6 ) 


(7) 

( 8 ) 

(9) 

( 10 ) 


All the model parameters re ported in Eqs. 1-10 are fitted, through a maximum likelihood 
unbinned algorithm (see § 3 in 1 Aiello et al.l 2014), to the Fermi -LAT (L 7 , z, T) data (see 


Tab. [[]). Among the three LF models, the LDDE model produces the largest log-likelihood, 
however a simple likelihood ratio test cannot be used to compare these non-nested models. 
We find that all three LFs provide an acceptable description of the LAT data (flux, luminosity 
and photon indices, see Fig. [[]), and more importantly predict comparable levels for the 
blazar integrated emission (see Tab. [[]). 


Blazars are known to have curved spectra when observed over a few decades in energy. 
It is thus important to have a reliable model of the high-energy component of the blazar 
SEDs. Here we use a double power-law model attenuated by the EBL: 


dN 7 

~dE 


= K 


EV (a / E x 7b ' 

E~J + \E b 


d —t(E,z) 


[ph cm 


-2r,-1 GeV -1 ] 


( 11 ) 


We rely on the EBL model of Finke et ah ( 20101 ). and use 7 a =1.7 and 76=2.6, which re¬ 
produces the long-term averaged spectra of bright BL Lacs with Ge V-TeV measurements 
(RBS 0413, Mrk 421, Mrk 501, see Aliu et al. 2012 : Abdo et al. 2011a b) and those of bright 
FSRQs (like 3C 454.3, 3C 279, 3C 273, etc.) observed by Fermi- LAT. Typically, all blazar 
spectra show a high-energy cut-off that reflects the distributions of the accelerated particles. 
This is located at E >lTeV and E <100GeV for BL Lacs and FSRQs respectively. Here, 
including such cut-offs makes very little difference because for BL Lacs the cut-offs are at 


2 Note that in Aiello et al. (12014 ) the exponents pi and p 2 were reported with the wrong sign. See Eq. |7| 


for the correct ones. 
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Fig. 1.— Observed redshift (upper left), photon index (upper right), 0.1-100GeV source- 
count (lower left), and 10-500 GeV source-count (lower right) distributions of Fermi -LAT 
blazars. For the upper panels, the continuo us so lid line is the PLE model convolved with 


the detection efficiency of Fermi -LAT (see Abdo et ah 2010a ). while for the lower ones it 


represents the predictions of the LF models. The 68 % uncertainty band in the lower right 
panel shows the prediction, for the 10-500 GeV source counts, of the LF and SED model. 
Error bars compatible with zero are 1 a upper limits for the case of observing zero events in 
a given bin. 


energies larger than those probed here, while for FSRQs, because of the larger average 
redshifts, the EBL efficiently suppresses their >50 GeV flux. 

For the model reported above, the high-energy peak is a function of Ef, alone, for fixed 
7 a and 7 &. We calibrated the relationships between Eb and the LAT-measured power-law 
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photon index via simulations and found that it can be approximated as logE{,(GeV) ~ 
9.25 — 4.11T (see left panel of Fig. [2]). The spectral curvature seen in bright LAT blazars 
is typically characterized using a logParabola model dN/dE oc /E / E 0 )~ a ~ l3iog ^ E ^ E °) (known 
to approximate blazar SED well only aroun d their peak), w here a is the photon index at 
energy Eq and [3 is the curvature parameter flNolan et al.l 2012). In order to ascertain that 
our SED model reproduces the correct amount of spectral curvature observed in blazars, 
we simulated LAT observations of ~1600 blazars with fluxes randomly extracted from the 
3LAC catalog and a spectrum described by Eq. [TTJ We treated these spectra as the real 
data and whenever the logParabola model was preferred over the power law at >4 a (as in 
Nolan et ah 20.12) we estimated the a and j3 parameters. As Fig. [2] (right panel) shows these 
are found to be in good agreement with the parameters of the real blazar set, validating our 
choice of the SED model. 

We thus use the above E b — T relation to predict the integrated emission of the blazar 
class that we compute as: 


r max—3.5 


x =6 


L™ ax =10 52 


EGB t- C/ 7 


(E 7 ) = 


dr 


d^ 


dL 7 - < h(L 7 , z, F)--^=r--j- [ph cm 2 s x sr 1 GeV 1 ], 


dE dz 


r min =i.o 


= in-3 

r. -J-VJ 


Lmin = 10 43 


( 12 ) 

where <I?(L 7 ,z, T) and are the LF and the spectrum reported above. Because the LF 
displays steep power laws at high redshift and luminosity, the only limit that matters is 
L™ m , which we set as the lowest observed luminosity. The normalization factor K of Eq. Qj] 
is chosen so that a source at redshift 0 and with index T, implying E h = E b (Y) given by the 
slope of Fig. [2] (left panel), has a rest-frame luminosity L 7 . We also make sure th at bo th the 


LF and SED models are able to reproduce the 10-500 GeV source counts (Ackermann et al. 


20131 ). which is important to obtain a robust estimate of the contribution of blazars to the 
high-energy EGB (see Fig. [[]). 

Integrating Eq. [12] above 0.1 GeV for the three LF models and averaging^ the results 
yields that all blazars (including the resolved ones) emit 5.70(±1.06) x 10 _6 ph cm -2 s -1 
sr^ 1 , where the error is dominated by the syste matic uncertainties (all similar in magnitude) 


on the Fermi -LAT detection efficiency flAbdo et alil2010 c). on the missing associations, the 
differences between the three LF models, and the scatter of the E b — T relation. When 
comparing this to the total EGB intensity of 11.3^} ® x 10 _6 ph cm -2 s~ x sr” 1 (AC14) we 
conclude that blazars produce 50^} 2 % of the total EGB. Since the resolved component of 
the EGB is 4.1(±0.4)xlO -6 ph cm -2 s _1 sr -1 (see AC14), and most of the detected sources 


3 We used a weighted average with l/erf (e.g. inverse of flux variance for each model) weights. 
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Fig. 2.— Left Panel: Simulated break energy Eb (for Eq. fill with 7 a =1.7, 7^=2.6) versus 
measured power-law photon index for a set of simulated blazars. The dashed line repre¬ 
sents the best fit described in the text. Right Panel: photon index (a, at 300 MeV) and 
curvature /3 (black data points) of the best-fitting logParabola models to simulated double 
power-law spectra (eg. Eq. [TU with 7 a =1.7 and 76=2.6). The gray datapoints show the 
parameters for all the bla zars (1 84) whose curvature is significantly detected in the 3LAC 
catalog (Ackermann et al. 2014dj). 


are blazars, we conclude that Fermi -LAT has already resolved ~70 % of the total blazar 
emission. 


3. Discussion 


Fig. [3] shows the spectrum of the integrated emission of blazarfl We find that the 
cut-off detected in the EGB spectrum is well explained by EBL absorption of the high- 
energy blazar emission. Above 100 GeV, the majority of the EGB is produced by nearby 
(z<0.5) low-luminosity hard-spectrum blazars. Below this energy, blazars cannot account 
for the entire EGB, in agreement with previous findings (see lAbd o et al. 20 10c) that, in the 
0.1-100 GeV band, unresolved blazars can account for only ~20% of the unresolved EGB 
intensity. Furthermore, it is difficult to accommodate a blazar population that produces a 
larger fraction of the <100 GeV EGB, because of the constraint placed by the level of the 


4 We neglected the secondary emission due to electromagnetic cascades created by electron-positron pairs 
generated in the interaction of y-rays with the EBL. 














































Table 1. Best-fit parameters of the LF models. Parameter values were computed as the median of all the best-fit 
parameters to the Monte Carlo sample, while the uncertainties (statistical only) represent the 68 % containment 
regions around the median values. Parameter names are reported as (k* or p*) and (£ or p£) for the PLE/PDE and 

LDDE models, respectively. 


Model 

Diffuse 3- 

A b 

71 

L* c 

72 

k * or pjj) 

T 

e or P2 

S 

z c 

a d 

M* 

p 

<7 

PDE 

k q«+0-19 

t>.»0_o 49 

-1 99+1*68 

2 - 80 ±o,25 


1 26"*" 0 ' 09 
i-ZD -0.08 

12 14+ 210 
1Z,A -1.74 

o 7q+0-56 
' y —1.30 

— 0 1 S+0-02 
U-1 -0.03 




9 99+O.O2 
^—0.02 

0 10+ 002 
U1 -0.02 

n 90+0.02 
°- 28 _0.01 

PLE 

K 7fi + °- 50 
'°-0.40 

i9-3i?;5 

3 i9+°- 51 
±y -0.40 

8 - 75 f2.42 

1 14+0-07 
— 0.08 

4 41+0-61 
— 0.64 

0 gi+O- 13 
uy± -0.15 

— 0 43+0 05 

o-^O-o 07 




9 99+O.O3 

-0.03 

0 10+ 002 
UiU -0.02 

O 90+O.O2 
-0.01 

LDDE 

5 4i+°- 57 

196+ 255 

iyD —130 

0 50~*~° -14 
U.OU-0 12 


-i oq+0.63 
1.83_0 35 

q qq+0.89 
3.3»_o 70 

3 16"*” 1 ' 45 

— 4 96"*" 2-25 
^• yo -4.76 

UD4 -1.05 

1 or + 0.19 
1ZO -0.17 

7 9 q + 2.17 
' z -2.99 

9 99+O.O3 
z z -0.02 

0 10+ 0-02 
ulu -0.02 

n 90+0.02 
u - 28 0.01 


integrated blazar emission (0.1—820 GeV), in units of 10 ® ph cm 2 s ^ sr , obtained by integrating the LF model between the limits reported in §0 
b In units of 10 —2 Gpc -3 . 

c In units of 10 48 erg s —1 for the PDE and LDDE models while units of 10 46 erg s —1 for the PLE model. 
d In units of 10 —2 . 
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small-scale anisotropies of the 7 -ray sky as measured by Fermi (Ackerman n et alj 2012). A 
blazar population that reproduces the 0.1-10 0 GeV source-coun t data flAbdo et al.l l2010ch 
can account for ~100% of the angular power (Cuoco et ah 2012 ), but only for ~20-30 % of 
the unresolved EGB. 


Therefore, the remaining ~50 % (~5.6xl0 -6 ph cm -2 s -1 sr -1 ) of the total EGB inten¬ 
sity (particularly at <100 GeV) must be produced by other populations or emission mecha¬ 
nisms. Star-forming galaxies and radio galaxies that, in addition to blazars and millisecond 


Jr -ml-u ^ -------- —j - -■ ^^ ? ------ ----- - --L—— . — —- - 

dio galaxies were shown to contribute 10-30 % of the EGB emission 

(Fields et ah 

2010 ; 

Makiva et ah 

2010 ; 

Ackermann et ah 

2012 

Inoue 

2011 ; 

Di Mauro et ah 

2013). By summing 

the contribution of star-forming galaxies h 

4ckermann et ah 

2012 ) and radio galaxies ( 

Inoue 


2 01 ll) to the contribution of blazars derived here (see Fig. |3]) we find that these three pop¬ 


ulations can naturally account for the intensity of the EGB across the 0.1-820 GeV range 
sampled by Fermi- LAT. This scenario does not ch ange if we adopt d i fferent models for 


the emission of star-f orming and radio galaxies (e.g. iMakiva et al.l 120101 : iFields et al.l 12010 


Di Mauro et ah 2013 ). 


This study shows that the source populations already detected by Fermi -LAT can ac¬ 
count for the entire measured EGB, leaving little room for other contributions. This can be 
used to constrain the emission from ‘yet undetected’ populations or diffuse processes. One of 
the most intriguing mechanisms that can produce a diffuse 7 -ray flux is the self-annihilation 
of DM present in the Universe. Indeed, if DM is composed of self-annihilating Weakly Inter¬ 
active Massi ve Particles (WIMPs) with masses of a few dozens to hundreds of GeV (see, e.g., 


Ber tone et ah 20051, for a review), then a diffuse GeV background may be expected from an¬ 


nihilations in DM halos across all cosmic epochs. This cosmological DM annihilation would 
thus contribute to the measured EGB, potentially yielding valuable information about the 
dark sector. No hints of a DM detection have been claimed up to now using the EGB. Yet, 


relying on the EGB intensit y (e.g. 


2014) or the anisotropy level (Gomez-Vargas et ah 


Abdo et ah 

20 

10 a; 

Bringmann et ah 

2014; 

Cholis et ah 


2014) 


Here, we use the main result of this analysis, that most of the EGB emission is produced 
by known source classes, to constrain the DM annihilation cross section. We rule out DM 
models that, together with point-like sources, overproduce the EGB emission at >2 a level. 


^Millisecond pulsars were shown to produce a negligible (<1 %) fraction of the EGB (e.g. Calore et al 


20 14) . 
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This is achieved by defining: 


X 


= Min 


' N 

E 


(F, 


EGB 


- sz/Fj 


i, ASTRO 


— Fi dm ) 2 (1 — £/) 2 


a. 2 


+ 


cr„ 


(13) 


where the sum runs over the N bins of the EGB spectrum. F iE QBi F^astro, F,dm are the 
intensities of the EGB, point-like sources and DM, srf is a renormalization constant of the 
nominal integrated source intensity and ov = {(Ji, astro / astro) its average uncertainty. 
In Eq. [13l ai is the sum (in quadrature) of the uncertainty on the unresolved EGB and the 
systematic uncertainty on the Galactic foreground (AC14). We use the uncertainties on the 
unresolved EGB because the uncertainties on the resolved source intensity are already taken 
into account in a^. The 2 a limits are found when the DM signal worsens the y 2 by >4 with 
respect to th e opt imized y 2 with a free DM signal normalization (and a free &/). Following 
Ackermami et al.l (12014a) . p redictions of the c osmological annihila tion signal were obtained 


using both the h alo m o del ([U lli o et al 


approach (jSerpico et al.l l2012t ISefusatti et al. 


2002 


2014). 


Fornasa et al.l 120131) and the power spectrum 


Though Eq. [13] neglects bin-to-bin 
correlations, we verified that our DM limits are within 10% of those obtained if, for each 
DM signal, we adopt the foreground model (from AC14) that gives the most conservative 
upper limit. 

An example of a ruled-out DM signal is reported in Fig. |31 while Fig. [I] shows the 
limits for DM annihilating to bb and t + t~ channels inclu ding their uncertainties due to 
the level of sub-halos in our Galaxy and in all DM halos (jSanchez-Conde fe Pradal I2014J; 
Ackermami et al.ll2014d). Our limits are co mpared to the conservative and sensitivity-reach 
limits reported in Ackermami et al. (2014c). The former assumes that the unresolved EGB 
is entirely due to DM annihilations, while the latter assumes the EGB is entirely produced 
by point-like sources and the DM annihilation is limited to what the uncertainties on the 
unresolved EGB allow. These represent extreme cases, while our limits represent a more 
realistic scenaric^, which, as expected, falls in between. Overall, our limits are 2-3 times 
more constraining than the conservative ones in A ckermann e t al. (2014c) thanks to the 
refined knowledge of the integrated emission from point-like sources derived here. 

This work shows that an analysis of the EGB and its components can constrain diffuse 
emission mechanisms like DM annihilation. The comparison of our limits with the sensitivity- 
reach limits of Ackermann et al. (2014c) shows that reducing the overall uncertainties is key 
to placing tighter constraints on DM annihilation. This can be achieved by refining the 
estimate of the emission from star-forming and radio galaxies. 


6 The sensitivity-reach limits in I Acke rmann et al. (2014cj) neglect the uncertainties in the integrated emis¬ 
sion of non-exotic sources which, once taken into account, will weaken the constraints on the cross section. 
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Fig. 3.— Top Panel: Integrated emission of blazars (with and without EBL absorption), 
compared to the intensity of the EGB (datapoints from AC14). Lower Panel: as above, 
but i ncluding also the emission from star-for ming galaxie s (gray band, I Ackerman n et ah 
2012 ) and radio galaxies (black striped band, Inoue 20111 as well as the sum of all non¬ 
exotic components (yellow band). An example of DM-induced 7 -ray signal ruled out by 
our analysis is shown by the solid pink line, and summed with the non-exotic components 
(long-dashed pink line). The inset shows the residual emission, computed as the ratio of the 
summed contribution to the EGB spectrum, as a function of energy as well as the uncertainty 
due to the foreground emission models (see AC14). 
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Fig. 4.— Upper limits on the self-annihilation cross section for the bb (top) and t + t 
(bottom) channels as derived in this work (see § [3]) 
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range of the theoretical predicted DM signal intensities, due to the uncertainties in the 
description of DM subhalos in our Galaxy as well as other extragalactic halos, adopting a 
cut-off mini mal halo mass of 1 CD 6 Mr.,. For comparison, limits reported in the lite rature are 
also shown (jAbramowski et al.ll201ll; lAckermami et abll2014al: lAleksic et al.ll2014l ). 









































